Analytical solutions for black-hole critical behaviour 
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Dynamical Einstein cluster is a spherical self-gravitating system of counterrotating particles, which 
may expand, oscillate and collapse. This system exhibits critical behaviour in its collapse at the 
threshold of black hole formation. It appears when the specific angular momentum of particles 
is tuned finely to the critical value. We find the unique exact self-similar solution at the thresh- 
old. This solution begins with a regular surface, involves timelike naked singularity formation and 
asymptotically approaches a static self-similar cluster. 
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General relativistic numerical simulation (numerical 
relativity) has revealed critical phenomena at the thresh- 
old of black hole formation in self-gravitating systems [H . 
When a parameter p, which parametrises a generic one- 
parameter family of initial data sets, is tuned to the criti- 
cal value p* , there appears a self-similar solution, which is 
called a critical solution. Beyond this value, the collapse 
ends in a black hole, its mass Mbh obeying the power 
law Mbh oc \p — p*\'^, where 7 is called a critical expo- 
nent. The critical behaviour is well described in terms 
of the behaviour of solutions around a self-similar solu- 
tion with a single unstable mode In this approach, 
self-similar solutions with regularity conditions are nu- 
merically found and a self-similar solution is numerically 
shown to be with a single, linearly unstable mode. See [3[ 
for a recent review of critical phenomena. Apparently, 
there still is a huge gap between numerical simulation 
and linear stability analysis. Moreover, one could sus- 
pect unresolved fine structure at the threshold because 
almost all results have been based on numerics with finite 
accuracy (cf. [1]). 

Here we show that there is a system where we can 
discuss critical phenomena in an analytical and exact 
manner. This is the spherical system of counterrotat- 
ing particles, first introduced by Einstein Q and later 
generalised to a dynamical case [1, Q ■ Using a coor- 
dinate r comoving to the radial motion of each shell, the 
line element in this spacetime is given by 

ds^ = -e^-^('^^'>dt^ + e^^^'^-^dr^-^R\t,r)dn^ (1) 

where dfl^ is the line element on the two-dimensional 
unit sphere. The Einstein equations and conservation 
law reduce 
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where ' = d/dt and ' = d/dr, M{r) and 2E{r) > —1 are 
arbitrary functions corresponding to the Misner-Sharp 
mass and the specific energy, respectively, and h{r, R) is 
given by 
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where L(r) is the specific angular momentum of counter- 
rotating particles. The energy density is 

M' 
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The motion of each shell is governed by Eq. 
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If we assume that the solution has a regular surface on 
which all regular metric functions and physical quantities 
are also analytic, this implies the Taylor-series expand- 
ability in terms of R^. The arbitrary functions M, E and 

then should be expanded as M = M^r^ -\- M^r^ + ■ ■ 
E = E2r^ + EiT^ + and = Llr^ + Llr^ • • •, if 
we choose r so that r = i? on the initial regular surface. 



R is expanded as i? = Ri{t)r -\- R^r^ 
surface, and using Eq. ([5]), we have 
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Observing the lowest order of Eq. ([3]) 
Ri{t) is given by 



• • • on a regular 
(8) 

the evolution of 
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If L| > 0, we can find that Ri{t) necessarily bounces 
and the neighbouring shell also does for sufficiently small 
r>0|i,i]. 

The case L4 = was studied in detail in Q. For 
this case, from Eqs. ^ and we find that an initially 
collapsing cloud inevitably form a central singularity af- 
ter a finite proper time. As for the motion of the shell 
for sufhciently small r > 0, it turns out that the quan- 
tity A = lim^^o -^('')/^^('') = Lq/M^ becomes impor- 
tant from Eqs. ^ and ([7]). Figure [T] shows the shape 
of potentials for different values of A. For A > 4, the 
region around r = necessarily bounces back, while, for 
< A < 4, the collapse continues to i? = 0. So, the end 
state of the collapse is a massless naked singularity or 
a black hole of finite mass depending on the value of A. 
The critical value A* is 4. The mass scaling law and the 
critical exponent will be described elsewhere [10|. 
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FIG. 1: The shape of the potential U for E = and L = AM 
with A = 3, 4 and 5. The horizontal axis is S = R/M. 

The self-similarity (homothety) requires that there is 
a vector field z such that C^gab = '^Qab- This condition 
implies for the line element ([1]) that there is a coordinate 
system (t, r ) such that i^, ip and S = R/r are functions of 
z = t/r [llj. For dynamical Einstein clusters, this means 
E{r) = E, M(r) = /ir and L{r) = Ir, where E, /x and I 
are constants. Then Eqs. ^ and ^ reduce the following 
closed autonomous system: 
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where a = v + hi\z\. We can eliminate z from the above 
equations and get 
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where the upper (lower) sign corresponds to the positive 
(negative) sign of dS/d{h\\z\) and 
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Equation pT|) reduces 
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where = e^^^^^dz^ = e^'^(dln |z|)^, ^ corresponding 
to the nondimensional proper time for the comoving ob- 
server. If we use as a coordinate in place of i, we get 
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where C = Inr and S plays a role of similarity variable. 
From Eq. (P, 
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If i?' = but R ^ 0, there appears shell-crossing sin- 
gularity and the model gets invalid beyond there. This 
condition reduces the following equation: 



STe''^y2{E-UiS))^0. 
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The solution is shell-crossing free where dS/d\n \ z\ < 0. 

When we rescale t and r as t = at and r — ar, we will 
get E = E, p, = /x/a, l = l/a, d- = a — Ina and S = S/a 
and the above equations are invariant. Hereafter, we fix 
this scaling so that fJ. = 1 and / = A, in which r coincides 
with the mass and is denoted as m. Note that we still 
have the freedom of rescaling as f = at because this only 
changes z to z — az and Eqs. (|10p - |fTT|) are invariant. Ex- 
cept for this gauge freedom, there are a three-parameter 
family of self-similar solutions parametrised by E, A and 
the initial value for (7 — fj j at S — Si. 

At the regular centre, m/R vanishes or S" = 00. To 
assure the circular constant being tt, ft,^/(l + 2E) ^ 1 as 
5 — > 00, implying E = 0. Moreover, we assume that the 
cluster has a critical angular momentum, i.e., the right- 
hand side of Eq. (fTT|) has a double root, implying A = 4. 
Then, we have specified two parameters and have a one- 
parameter family of self-similar solutions. For E — 
and A = 4, the shape of the potential is shown in Fig. [TJ 
Equation p4)) then has essentially three solutions and 
all others are generated through time reversal or time 
translation of those solutions. The first one is a static 
solution at the top of the potential. For this case, we 
assume ip and R have no dependence on t but do not for 
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V. This is because djdt may not coincide with a timeHke 
KiUing vector. Then, Eqs. pO)) " ((TT|l yield the following 
solution: e^" ~ Cglzj, e^"^ — 32, and 5* = 4, where cq is 
a constant and set to be unity, using the rescaling of t. 
The resulting metric is given by 

ds2 ^ -RdT"^ + 2(ii?2 + B^dV?, 

where we have implemented the coordinate transforma- 
tion: dt^l\t\ = MT^ and 4to = R. The density e is 
given by 47re = l/(64m2) = l/(4i?2), This spacetime 
has a timelike naked singularity at the centre and suffers 
from a solid angle deficit so that the area of the sphere 
R = const divided by the squared radius is not An but 
27r. It can be easily shown that static singular solutions, 
which were discovered in [l^ , all fall into this self-similar 
static solution after an appropriate coordinate transfor- 
mation. 

The remaining two solutions are dynamical. The sec- 
ond one is the solution, in which S begins with infinity at 
^ = — oo, monotonically decreases and asymptotes 4 as ^ 
increases to oo. The third one is the solution, in which 
S begins with 4 at ^ = — oo, monotonically decreases to 
at a finite value of For these two solutions, we can 
express e'^ in terms of S using only elementary functions 
as an integral of Eq. p2|) : 
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where A is an arbitrary constant. This is actually a sub- 
class of solutions whose expressions were given in terms 
of elementary functions in [9']. Only the second solution 
has a regular centre because m/ R is always equal to and 
larger than a quarter for the static solution and the third 
solution, respectively. 

Let us concentrate on the second solution. Around the 
regular centre S — oo, from Eq. (HH]), we have 

^^S^'^ + A^S^'^ + A + 0{S-^'^). (19) 

Then, from Eqs. (jlOp and (fTTjl . we have for the lowest 
order term S « (9/2)i/3x2/3|z|2/3 and « X\z\, as 
1^1 — > cxD, where the constant X can be set to be unity by 
rescaling t. This is exactly the behaviour of the regular 
centre: R cx m^^^ ^nd e"^ ^ 1 at m — > 0. In this regime, 
^ — z-l- const, and t coincides with the proper time at the 
regular centre. Equation implies 47reo(t)t2 = 2/3, 
where eo(t) is the density at the regular centre. The 
central density diverges to infinity at f = 0, resulting in a 
central singularity. In other words, i = is characterised 
with the appearance of the central singularity. Let us 
assume the expandability for the density on the regular 
surface t = const < 0: 



(20) 



On this surface, we can write down S = R/m as 

1 / ^, \ -1 
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Substituting Eqs. ^ and into Eq. and com- 
paring both sides, we get ei = 0, €2 = (160/3)7re2 and 
£3 = 16-y2/37r3/2^£^/2_ Therefore, the density around 
the regular centre can be expanded as 

1 40 i?2 4 R3 
e= T + rH A— + 0{R^). (22) 

The density takes a local minimum at the centre. If we 
assume analyticity on the regular surface, A = is con- 
cluded and there no longer appear odd power terms of R 
in the expansion. This analyticity requirement has been 
imposed for the identification of critical solutions 0]- So 
we will identify the second solution for A = with the 
unique threshold solution. Note that the analysis below 
is nevertheless also applicable to the case A ^ 0. 

Let us see other physical properties of this solution. We 
can find from Eq. that e'^ can be zero. Let a{> 4) 
such that e"' = at 5 = a. a(= 4.8133 • • •) is a root of 
the following transcendental equation: 



VS{S + 12) ^ 12 In 



VS + 2 
VS-2' 



Around there, from Eqs. (jlip and ([H]), we can find \S — 
a\ cx |zp and e°' oc \z\^ as z — > 0, where f3 = a2/(a2-|-16) 
and hence 1/2 < /3 < 1. In this regime, ^ = =fB|z|'' -|- 
const, where i? is a constant. This means that it takes 
only a finite proper time to reach S = a or t = 0. Thus, 
t = is only a coordinate singularity and the solution can 
be extended regularly and uniquely beyond t = 0. This is 
clear in Eq. (|15|). where there is no singularity at 5 = a. 
The density profile at t = follows an exact power law 
47re = l/{aR'^). Because of self-similarity, 2 = also 
corresponds to infinity (|t| < 00 and m = 00). This 
really corresponds to the surface of infinite area R — 00. 
It follows that the density falls off as Aire « l/(ai?2) as 
i? — > 00. This behaviour is time- independent. 

Each shell dynamically approaches the specific radius 
5 = 4 as ^ ^ 00. For 5 — > 4, from Eq. pS)) . we have 
« -41n(S' - 4). Then, from Eqs. ^ and HI]), we 
have {S -A) w g-'^l^l'^" and = AC\z\^/^ as \z\ 
00, where C > is a constant. This means that the 
asymptotic solution is the static solution. In this regime, 
^ « 8B|z|i/2 _^ const, and r w -8mln(S' - 4). Hence the 
collapse approaches the static solution with an infinite 
proper time. 

Since dS'/(iln|z| > for t < for the threshold solu- 
tion, we need to check whether there is a root in (a, 00) 
of Eq. (|17p or the following equation: 
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When S increases from a to cx), the left-hand side in- 
creases from 0, takes a maximum and decreases to 2/3. 
This maximum value is 0.76412 • • •, well below unity. 
Since dS/d\n\z\ < for t > 0, the threshold solution 
is free of shell-crossing. 

Figure [5] shows the evolution of the density in terms 
of S from Eq. p6)) . As time proceeds from < = — oo 
to t = cxD, S monotonically decreases from oo to 4 and 
the density e observed by a comoving observer monoton- 
ically increases from to the value for the static solu- 
tion l/(2567rrn^). It smoothly crosses z = a,t S — a 
and e = l/(47ra''m^). Figure [3] shows the ratio of the 
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FIG. 2: The nondimensional density A-Kem in terms of S for 
the exact threshold solution. 

local density e to the averaged density inside the shell 
m/(47ri?^/3) , in which the horizontal axis is I /Vs. This 
ratio must be unity at the regular centre and 1/VS is 
proportional to the area radius around the regular cen- 
tre. We can see that the density takes a local minimum at 
the regular centre S' = ooor2; = — oo, increases around 
the centre, takes a maximum and decreases outside. At 
spatial infinity S = a or z = 0, this ratio becomes 1/3, 
implying that e falls off in proportion to R~^. The ra- 
tio goes below 1/3 as z increases further and comes back 
to 1/3 at S* = 4 or 2; = oo, implying that e again gets 
proportional to as t —> oo. 

There is no trapped surface because S — R/m > A is 
satisfied everywhere. So, the central singularity is naked. 
In fact, this spacetime is a member of solutions for which 
the causal structure is shown in Figure 1 of [l3|. The 
central singularity is naked and timelike. 

The general solutions of the dynamical Einstein clus- 
ter is exactly solved in terms of elliptic integrals using 
the mass-area coordinates [13] and here the critical 
self-similar solution is uniquely obtained in terms of el- 



ementary functions. This system provides a tractable 
laboratory for studying how generic critical collapse ap- 
proaches the threshold solution in both linear and non- 
linear regimes. 
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FIG. 3: (47reJ?^)/(3m), the ratio of the local density to the 
averaged density interior to the shell in terms of l/\fS for the 
exact threshold solution. 
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